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ABSTRACT 

We report on a novel method to solve the basket- weaving problem. Basket- weaving is a technique that is used to remove scan-line 
patterns from single-dish radio maps. The new approach applies linear least- squares and works on gridded maps from arbitrarily 
sampled data, which greatly improves computational efficiency and robustness. It also allows masking of bad data, which is useful for 
cases where radio frequency interference is present in the data. We evaluate the algorithms using simulations and real data obtained 
with the Eff'elsberg 100-m telescope. 
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In single-dish radio astronomy receiving systems often provide 
only one or few pixels because the necessary feed horns cannot 
be as densely packed as for example pixels on a CCD in op- 
tical astronomy. Therefore, to obtain sufficient mapping speed, 
most observers choose to observe in so-called on-the-fly mode, 
also known as drift scanning. Here, the telescope is continuously 
driven across the target area, forming scan lines or even just a 
single continuous path (e.g., following a Lissajous curve; see 
alsol^vacs 2008). This minimizes the duty cycle and also al- 
lows very shallow observations by scanning the map as fast as 
necessary, constrained only by the maximal telescope speed. 

Because gain of the receiving system, the atmosphere, and 
contribution of the ground may vary with time, one often en- 
counters image artifacts in the form of stripes (the scan pattern) 
after regularly sampling the observations on a pixel grid (grid- 
ding). One poss ible so l ution to remove such effects was pro- 
posed by Hasla met al.l (Il970h . called the basket- weaving tech- 
nique, as an "iterative process which minimized the temperature 
differences at crossin g points". Subsequen tly, this was applied to 
the 408-MHz survey (iHaslam et al.lll98lh . 

Basket- weaving refers to a strategy where the target area of 
interest is mapped twice with (approximately) orthogonal scan- 
ning direction. Since the observing system should receive the 
same contribution from the astronomical sky, the difference of 
both measurements should only be due to the same effects caus- 
ing the image artifacts mentioned above. The aim is to find the 
best correction (offsets and/or gain factors) to each scan line un- 
til both maps lead to a consi stent result. 

For the 408-MHz survey 'Hasl am et all (11981 ) optimized for 
moderate-gain variations of the receiver using a linear function 
for each scan line. The optimization itself was applied to the 
scan lines one by one and the procedure was repeated until the 
residual effects were of the same order as the noise level. 

There also exist other methods to deal with scanning effects 
that do not rely on basket- weaving. For example, Sofue & Reich 
(1 19791) described a filtering approach (based on unsharp- 
masking) to suppress artifacts that usually have a different spa- 



tial scale — relatively elongated in one dimension and smaller 
than the beam size in the other (at least if the maps were spa- 
tia lly fully- sampled). Anothe r filtering technique was proposed 
by lEmerson & Graev3 (Il988h . where again two orthogonal maps 
are necessary that are combined in the Fourier space using a 
weighted averaging. Both methods have in common that the spa- 
tial scales suppressed by the filter have to be known or approx- 
imated beforehand. Removing certain scales may affect parts of 
the astronomical signal as well. We point out that these filters can 
only remove additive effects, e.g. emission of the atmosphere, 
and not multiplicative ones, e.g. caused by gain variations or at- 
mospheric dampening. 

More recent examples of large projects that employ ba sket- 
weavin g can be found in Wolleben et al. (2010) and Peek e t al.l 
(1201 lb . Since the data volume produced by modern receiving 
systems has substantially increased, it is often computationally 
inefficient to follow an iterative approach. This is especially the 
case for large-area surveys carried out with multi-feed systems, 
like the Effelsberg-Bonn H i survey (EBHIS, Winkel et al. 2dTc| 
lKerpetal.ll2Qnl) . where about 320 individual scan lines (per 
coverage of a 5 x 5 deg^ field) are produced. For each scan line 
and iteration the residual RMS had to be computed for ~ 10^ 
pixels in each of thousands of channel maps. 

In this paper we present an algorithm that is based on linear 
least- squares for solving the basket- weaving problem directly. 
Furthermore, the problem matrix has to be computed only once 
per map, allowing one to re-use it for each channel map in the 
case of spectral data sets. Another clear advantage is that the al- 
gorithm operates on gridded data, such that no interpolation is 
need ed to find the inte nsities on the intersection points (com- 
pare Haslam et al. 1981). The algorithm still is only able to re- 
move additive effects (due to the linear least- squares approach) 
which, in practice, rarely imposes a limitation as modern re- 
ceivers are relatively stable such that gain variations are usually 
small. Furthermore, using calibration techniques (e.g., keeping 
track of the signal of a noise diode) gain effects imposed by 
the receiver can easily be handled. Another phenomenon with 
both a multiplicative (attenuation of the astronomical signal) and 
additive effect (emission of the atmosphere itself) on the astro- 
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nomi cal signal is opac ity, r, of the atmosphere. Using models 
(e.g. JPardo et alJ '2006') and simultaneous measurements of wa- 
ter vapor in the at mosphere (e.g. using a water vapor radiometer, 
iRottmann & Roy|[200 7) one is able to determine t(0 with suffi- 
cient accuracy to remove most of the associated gain and offsets 
effects. However, residual uncertainties will manifest themselves 
mainly as baseline offsets and can be dealt with using our ap- 
proach. 

Techniques v ery similar to bas ket-weaving exist. For the 
Planck mission (iTauber e t al.' '201 Q[), for example, a me thod 
called de-striping (e.g., lDelabrouillll998|; lEfstathioul[2QQ7l and 
references therein) was used to remove the differences between 
the individual stripes of data obtained by observations during 
a slow precession of the Planck satellite. This leads to a scan 
geometry quite different from typical basket- weaving scan pat- 
terns. As a consequence, one "scan line" can have intersection 
points with all others, which increases the complexity of the 
problem even further. 

The paper is organized as follows. In Section[2]we show how 
to solve the basket- weaving problem using a linear least- squares 
approach in general. As previously mentioned, one can find an 
implementation that operates on gridded data, which can greatly 
improve computational efficiency if the data were spatially over- 
sampled (i.e., the number of intersection points is much larger 
than number of pixels in the map) or if spectral data with many 
channel maps is to be processed. This is presented in Section[5] 
and accompanied by simulations. Two examples using real data 
are discussed in Section|4]and we conclude with a brief summary 
(Section©. 

2. Least-squares basket-weaving 

As discussed in the introduction, basket-weaving aims to mini- 
mize the difference, dp, of the measured intensities, y^J^^ and yj^, 
at each intersection point (j, /) in two quasi-orthogonal cover- 
ages. Here, / and j are indices over the number of scan lines, /^^^ 
and J^'^\ in the first and second coverage. In the following we 
will omit the superscripts (1) and (2) when it is clear from the 
context which coverage is meant. In Fig.[T]the geometry of the 
basket- weaving problem is schematically shown. It is not nec- 
essary that the scan lines are plane parallel, as long as the two 
coverages have a sufficient number of intersection points. One 
could even use very complex scan patterns, for example based 
on Lissajous curves (see also Kovacs 2008) . 

Since the astronomical emission, Tjj, should be the same in 
all coverages, we find 



= pr-P 



(2) 



1.../, J = 1.../. 



(1) 



Thep^^^ are the offset values for coverage k and the ^-th scan 
lin^l The latter are unknown quantities, we only know their pair- 
wise difference. The aim is to determine a set of p^^^ that simul- 
taneously fulfills all / X / equations. To enable a least- squares 
approach, it is useful to find a matrix notation, 

^ ^VrlJ\,{r mod J) ^VrUHr mod /) 



^ {^s,{r mod /) " ^s,VrlJ\+l)Ps 
s=l...I+J 



(2) 



s=l...I+J 



^ In Section [331 we also discuss the more general case of arbitrary 
(smooth) offsets which can be parametrized by polynomials. 
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Fig. 1. Geometry of the basket-weaving problem. Usually two 
(quasi) orthogonal maps are observed in so-called on-the-fly- 
mode (OTF) where the telescope constantly drives across the 
area. Each of these scan lines is visualized by an arrow (red: 
first coverage in longitudinal direction, blue: second lateral cov- 
erage). At the crossing points dj^i the contribution of the true 
brightness temperature should be the same for both coverages, 
while baseline (offsets) may differ. The distance between two 
recorded samples (dumps) on a scan line is often smaller than the 
separation of the scan lines. Depending on the projection and co- 
ordinate system chosen, the scan lines may not be plane-parallel. 



(1) (2) (2) 



pf) and 



with = (p{,p!^) 

= {d[, d^,. . .,dj) and dj = (dij, d2j, . . . , djj). The notion 
[xj is used for the floor of x (the largest integer not greater than 
x). The matrix A has the size (/ x /, / -h /). In short, 



D=AP. 



(3) 



This linear equation can be solved using linear least- squares 
(e.g., via singular value decomposition, SVD) by minimizing 



\AP - D[ 



(4) 



Different noise realizations usually lead to different least- squares 
solutions. Furthermore, the above equation is degenerated (i.e., 
has no unique solution intrinsically). 



2.1. Degeneracy 



For example, one could just add an arbitrary value p to all p] 
and p^P without changing dp. Therefore, it is necessary to ap- 
ply additional constraints. One common way to do that is to in- 
clude dampening (of the parameter vector) by minimizing the 
weighted sum of squared norms instead 



\AP- 



Df -h \Pf , 



which is equivalent to 



(5) 



(6) 



with the dampening v alue A. This will lead to a bounded solu- 
tion P (see for example lBovd & Vandenberghel2004l:lBussll2009L 
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and references therein). The choice of a reasonable value for A 
depends on the problem itself. We will discuss this briefly in 
Section[3?6l 



3. Improving computational efficiency 

The method presented in Section[2] is computationally demand- 
ing, because each intersection point requires one additional row 
in the linear equation. Furthermore, one has to determine the in- 
tersection coordinates, (a(j, /), S(j, /)), and potentially interpolate 
the data to obtain the values yjj which usually will not match the 
spatial position, (of^, Sa), of the recorded samples. 

Often the recorded data are more densely sampled than re- 
quired by the Nyquist sampling theorem. For these cases we 
propose a modified basket-weaving algorithm that operates on 
gridded data. We will also show that this approach has several 
additional advantages (and few disadvantages). As a prerequi- 
site, a convolution-based gridding algorithm is briefly discussed, 
which can be implemented to allow serial data processing (use- 
ful for large data sets) and can be easily parallelized. 

The de-striping algorithm used for the Planck mission is very 
similar to the least-squares technique discussed in Section^ It 
simplifies the problem of interpolating the raw data onto the 
intersection points by working on a special pixelization of the 
sky — the HEALpix grid ( Gorski et al. 2005) — which is well- 
suited for the special scan geometry. 
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Fig. 2. On top of the geometry sketched in Fig.[T]a regular pixel 
grid (black lines) is constructed. In this example there are fewer 
pixels than intersection points. The pixel grid does not need to 
be aligned with the underlying scan lines. 



Eq.O 



3.1. A convolution-based gridding method 

To grid a number of data points ja^cta^ ^a) into a map consisting 
of M X N pixels {m,n) with the world coordinates (of^,«, ^m,n), 
one can apply the formula 

RmAy^ = 7 7 ' w) 

where a is an index over the total number of measured data 
points to be gridded and w{am,n ~ cta^^m,n - ^a) is a weight- 
ing function (or so-called convolution kernel). For the remain- 
ing part of this paper we will use a Gaussian convolution kernel 
that has convenient Fourier properties but slightly degrades the 
spatial resolution in the resulting mapE|. Now (m, n) is not the 
crossing point of scan line / (first coverage) and j (second cover- 
age), but refers to the pixel with (pixel) coordinates (m, n) in the 
gridded map. 

In the following we make use of the shorter notation 

W(m, n\ a) = W(am,n - OLa. ^m,n " ^a) and Wm,n = Z« ^'^ 



3.2. Re-sampling the problem matrix 

If the two coverages were (independently) gridded to yield the 



maps and we can compute a diff'erence map (see also 
Fig. [2] for a sketch of the geometry), similar to the approach in 



r>(2) 



^ Note that the size of the kernel is loosely linked to the beam size of 
the instrument. The resolution in the gridded map must be good enough 
to ensure Nyquist sampling, but it should also be fine enough such that 
the gridding kernel is fully sampled. Therefore, to avoid extreme over- 
sampling of the image, the gridding kernel should not be much smaller 
than the instrumental resolution. In this paper we work with a kernel 
having half the beam size, leading to a degradation of about 12% of the 
spatial resolution in the final map. 



= RmAT^'^] + RmAP^'^] - RmAT^^^] " ^mAP^^^] (9) 

(10) 



^m,n ^(1) 



where we have used that = which one can 

safely assume if both maps are fully sampled and the true sky 
brightness has not changed. 

As mentioned, a is the total number of measured data points 
being gridded (per coverage). For later use, it obviously makes 
sense to also store the scan-line numbers / (coverage 1) and j 
(coverage 2) each dump a belongs to, as well as the number of a 
dump in its scan line. The latter shall be denoted with u (cover- 
age 1) and V (coverage 2). Since the number of dumps per scan 
line may vary, we call Ui iVj) the total number of dumps in scan 
line / (j) in coverage 1 (2). 

When we again assume that each dump in a scan line is sub- 
ject to the same off'set (i.e., we work with one off'set per scan 
line), we can write 



dm,n = + —(^ ^ pf'' ^ w ((m, n\ a^^\i, u)) 



^m,n i=l u=l 
V 



^m,n j^l v=l 



(11) 



As before, we obtain a matrix notation by just flattening dm,n 
using the mapping (m, n) ^ r with m{r) = (r mod M) and 
n(r) = [r/MJ. As before, the matrix is still only dependent on the 
geometry of the two coverages (i.e., the sky positions at which 
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50 100 150 50 100 150 

Column Column 



Fig. 3. Left panel: Basket- weaving matrix A for a small test case (2 x 100 strictly orthogonal scan lines, gridded into a map of 
15 X 15 pixels). The left/right part of the matrix (columns 1 to 100/columns 101 to 200) links the offsets (100 each) associated with 
first/second coverage to the difference map (flattened to form the vector D). Right panel: For the same example the coordinate grid 
was rotated by 45°. Evidently, now different offsets contribute to a pixel compared to the left panel. 



data were taken) and the choice of a gridding kernel. 



S^l...I+J 



^^^^ 

m{r),n{r) 



— ^r,sPs 



s=l...I+J 



with 







(2) 
m{r),n{r) 



a < s < , 
else 



(12) 



(13) 



and i(s) and j(s) being the mapping functions to obtain from 
the sequence s the appropriate scan line number for the two 
coverages, i.e., i(s) = ^ for 1 < 5- < / and j(s) = s - I for 
I-\-l<s<I-\-J. Again, Ps is the vector containing all offsets 
as defined in Section[2j 

The two necessary quantities, D, and the matrix, A can be 
easily obtained using the following recipe 

1. Grid all data points of coverage 1 and 2 into two separate 
maps, 7?^,'^^ and keep the weight maps w^'^\ 

2. Subtract both maps according to Eq. dS]) and flatten the result 
to obtain D. 

3. During the gridding process, compute additional / + / maps: 
for each of the / (/) scan lines in coverage 1 (2), set all data 
points to 1 and grid separately each scan line into one indi- 
vidual map. Divide the first / maps by w^^l^^n ^nd the following 
/maps by w^,^„. 

4. Flatten each of the resulting I -\- J weight maps; the resulting 
vectors form the columns of the matrix A. 

While equation Eq. ([TT1) looks more complicated than 
Eq. ([B, the problem matrix is just a "smoothed" version of the 
former, having the same basic structure. Likewise, the calcu- 
lated offset values Pg are just a "smooth" estimate of the true 
offset values, having the same effect on the gridded data. The 



"re- sampled" matrix A has MxN rows and / + / columns, i.e., it 
can be of much smaller size if the the original maps were over- 
sampled. 

In Fig. [3] we have visualized the basket- weaving matrix for a 
small example. Each row in the matrix corresponds to one pixel 
in the difference map, i.e. it links various scan line offsets (hav- 
ing different weights) to it. The left part of A refers to the first 
coverage, while the right part refers to the second (here values 
in A are negative to achieve subtraction of both coverages). The 
example in the left panel of Fig. [3] is for strictly orthogonal scan 
lines. To demonstrate a more complicated case, we simply ro- 
tated the coordinate grid by 45° (but not the scan lines), leading 
to the matrix in the right panel of Fig.[3l 

As before, additional constraints are necessary (e.g., the 
dampening strategy), as the problem is still degenerated. 

3.3. Calculating the correction map 

Using the gridding-based approach facilitates calculating a cor- 
rection map without the need to re-run the gridding process (with 
the computed offsets applied). This is possible because the off- 
sets are additive, and the matrix A already contains all necessary 
weight factors. We find that 



^s=l m{r),n{r) ''"^ ^s-l m{r),n{ry ''"^i 



m{r),n{r) m{r),n{r) 



(14) 



m{r),n{r) m{r),n{r) 



is the overall correction map (after re- shaping to two dimen- 
sionsfl If one similarly applies \A\ to subsets of P, one can also 
obtain correction maps for the individual coverages. 

This makes the proposed method extremely useful when 
a large set of observations with identical spatial positions of 
the data samples have to be processed, as for example in the 
case of spectroscopic observations or state-of-the-art continuum 



^ Taking the absolute value of the elements is necessary since the 
matrix A contains negative values in the columns / + 1 to / + /. 
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Fig. 4. Simulated maps. The left panel contains the input sky model. From this, two coverages were drawn and noise was added as 
well as a random offset for each scan line. The middle panels show the average of both maps (for the clean and dirty case, the former 
with and the latter without offsets applied), the right panel the difference dm,n according to Eq. dS]). Note that different intensity scales 
were used. 
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Fig. 5. Solving Eq. ([T2l) for the difference map dm,n (Fig.|4]right panel), one obtains a solution for the offset values P. For comparison 
we show in the left panel the reconstructed difference map (which is basically a noise-free version of the input difference map). Via 
Eq. ([T4l) one can easily compute a correction map (second panel) that can be subtracted from the original "dirty" map (see Fig. [4] 
third panel). The result is a cleaned data set (third panel), from which image artifacts are visually removed. The right panel shows the 
difference of the clean and cleaned map. Residual stripes appear, but they are suppressed by about an order of magnitude compared 
to the correction map (note the different intensity scale). 



measurements (with many channels). For these it is likely that 
the offsets P will be a function of frequency, such that basket- 
weaving has to be applied to many spectral channel maps. Then 
the matrix A is the same for all channel maps, only D needs to 
be computed for each channel, which is usually done in the grid- 
ding process anyway. 

3.4. Simulations 

We performed simulations to test the proposed basket-weaving 
method that also allowed us to analyze the quality of the solution 
and the performance. The simulated map parameters resemble 
those of the Effelsberg-Bonn HI Survey (EBHIS, Winkel et al. 
l2QlQl:lKerpetal.ll201lh since this is one of the potential appli- 
cations of least- squares basket- weaving. For EBHIS the angular 
resolution is 10', which means that for a typical field size of 
5° X 5° one needs about 100 x 100 pixels of size 3' to achieve a 
fully sampled map after gridding. The number of observed scan 
lines is higher with about 320 scan lines per coverage, which 
means that the number of intersection points is an order of mag- 



nitude larger than the number of pixels in the final map. Each 
scan line contains about 160 data points (dumps). 

For our simulations we used a simple sky model containing 
several Gaussians on top of a (2-D) polynomial background; see 
Fig. m (left panel). From this sky model we sampled two maps, 
one with longitudinal and one with lateral scan direction. To 
both of these coverages we added Gaussian noise and an offset 
to each scan-line, the amplitudes of which were sampled from 
a Gaussian distribution (with similar standard deviation as the 
noise distribution). Figured shows the average (middle panels) 
and difference (right panel) maps, clearly containing a strong 
pattern caused by the offsets. 

After constructing and solving the matrix equation Eq. ([T2b 
based on the difference map in Fig. |4l (right panel), one obtains 
a least- squares solution for the offsets P. Multiplying the ma- 
trix A with P leads to the reconstructed difference map (Fig.O 
left panel), which is basically a noise-free version of the for- 
mer. Likewise, with Eq. ([T4l) we can compute a correction map 
(Fig. [3 second panel) that leads to a cleaned average of both cov- 
erages (Fig. [51 third panel) when subtracting from the dirty av- 
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erage map (Fig.lH third panel). In this cleaned map, visually no 
image artifacts (stripes) appear. However, the difference between 
the clean and cleaned data shows some residual errors (Fig.lH 
right panel), which are about an order of magnitude smaller than 
the original image artifacts/stripes. 

3.5. Polynomial offsets 

Up to now we only accounted for the relatively simple case of 
constant offsets for each scan line. Of course, in reality it might 
well be that there is some underlying trend in the data of a 
scan line, e.g. due to weather effects or variable radiation from 
the ground. Then one possibility that still allows a linear least- 
squares approach is to parametrize the drifts using polynomials. 

Usually, individual data samples in a scan line will have the 
same integration time, such that one could construct the polyno- 
mials as a function of dump number per scan line. To obtain a 
similar order of magnitude for all polynomial coefficients (up to 
the maximum order of Op), it is possible to re-scale the integer 
dump number I <u <Ui (and 1 < v < Vj) to the interval [0, 1] 
by dividing the dump number by the total number of dumps in a 
scan line. Eq. ([T]) then reads as 

<'«=Zrf,;(ff-E;'3fef. (15) 

However, with this approach one has to parameterize and in- 
terpolate the position of each data point on the scan lines to ob- 
tain Ui and Vj, which usually are not integer numbers. This adds 
complexity to the problem. With the gridding-based method, this 
will naturally be solved, since the gridding is equivalent to inter- 
polation. 

One can easily convert Eq. ([TSl) into a matrix equation similar 
to Eq. (|2l). The matrix A will then not only contain entries + 1 and 
-1 but also the values of (ui/Uiy and (Vj/Vjy, respectively. It 

has / X / rows and (Op'' -h 1) x / -h (of'' + 1) x / columns. 

Likewise, the gridding-based approach can be adapted to al- 
low for higher polynomial orders. Here, we omit the (lengthy) 
equations, but just extend the recipe given in Section [T2l Only 
step (3) needs to be slightly changed 

3. In the gridding process, compute additional (of^ + I) x I -\- 

(of^ + 1) X / maps: for each of the / (/) scan lines in cov- 
erage 1 (2), set the data points to (u/Uiy and (v/VjY, re- 
spectivel}0, and grid these data separately for each scan line 
and polynomial coefficient into one individual map. Divide 
the (of^ -h 1) X / maps associated with the first coverage by 
(of^ + 1 ) wL\^„ and the other (of + 1 ) x / maps by (of + 1 ) w^,^ . 

Again, the resulting matrix will have the same number of 
columns as for the non-gridding-based approach, i.e., (of -h 1) x 

/ + (of + 1) X /, but a different number of rows, M xN. 

Note that also the reconstruction formula must be slightly 
changed to account for the increased number of parameters that 
contribute to each pixel in the map. This is done by simply 
changing the denominator of Eq. ([T4l) to (of + l)wJ;J^^„ + (of + 

l)Wm,«. 

Of course, one cannot solve the system for an arbitrarily high 
polynomial degree. Artificially increasing M and N (the spatial 

^ Here u and v are again the integer dump numbers and not interpo- 
lated quantities. 



resolution of the gridded map) will not help either, because ad- 
jacent pixels will be correlated and not add information. In total 
M X N should be much larger than I -\- J > M -\- N, which just 
means that a sufficient size of the map is necessar}0. 

One could also choose a completely different set of basis 
functions to parametrize the offsets, not only simple polynomi- 
als. Furthermore, it is possible to choose different free parame- 
ter(s), e.g., elevation instead of dump number. Both are possible 
because one just needs to grid the values of the functions for each 
coefficient into the matrix — the structure of the basket- weaving 
equation is still the same. 

In Fig. [6] we show results for higher-order polynomial off- 
sets. In reality it is not necessarily clear what the best recon- 
struction polynomial degree would be. Therefore, not only the 
cases where the simulated polynomial offsets match the poly- 
nomial degree used for the basket- weaving matrix are shown. 
For example, the first column shows the results for using vari- 
ous reconstruction degrees, all applied to data containing con- 
stant offsets for each scan line. Due to the increasing number of 
parameters — using a third-order reconstruction polynomial re- 
quires already ~ 2500 parameters, which is comparable to the 
total number of pixels (10"^) — the solution becomes increas- 
ingly worse. The same is true if higher-degree polynomials are 
used to generate the data (from left to right in each row), but 
the reconstruction polynomial order is not sufficient to describe 
these. But also along the diagonal, where both polynomial de- 
grees match, the results become worse for increasing degrees. 
Again, this is a question of number of free parameters, which 
are less constrained to the lower right (the noise level was kept 
constant for all cases). 

Nevertheless, in all cases the results are still better than with- 
out basket- weaving (the top row of Fig. [6] contains the dirty im- 
ages). 

3.6. Choice of the dampening parameter A 

In Section im we discussed how the (degenerated) basket- 
weaving problem can be solved by applying a dampening 
scheme that constrains the solution. However, it is not directly 
evident how large A should be. For our example, using first-order 
polynomials, we processed the data with a wide range of differ- 
ent A values. To improve statistics, this was repeated 30 times, 
each time using different noise realization and offset coefficients. 
The input model was subtracted from each reconstructed map to 
calculate the residual RMS, which we plot in Fig. [7] (normalized 
to the RMS in the clean data set residuals). For reference, the 
plot also contains a horizontal line (green) for the value of One 
and the normalized residual RMS for the dirty maps. The for- 
mer defines the noise level for perfect reconstruction, the latter 
resembles the case where no improvement was achieved com- 
pared to the dirty maps. 

Fig. [7] reveals that a wide range of A values (three orders of 
magnitude) leads to a consistent solution. For large A, which 
means that the constraint on the offsets is tighter, the resid- 
ual RMS converges to the value of the dirty images. Here, the 
solution is in fact bound to zero. Toward small A the residual 
RMS constantly increases, a sign for divergence of the solution 
because the degeneracy of the problem is not sufficiently sup- 
pressed by the dampening anymore. 



^ The equality is valid if the scan lines were perfectly Nyquist sam- 
pled. In all other (oversampled) cases, the grid-based approach will re- 
duce the size of the matrix and increase computational efficiency. 
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Fig. 6. Result of the least- squares basket- weaving correction for different polynomial degrees. The four columns belong to the 
generated polynomial offsets (from constant offsets in the left column to a third-order polynomial in the right column). The top row 
shows the dirty maps for these. The four bottom rows refer to different reconstruction degrees, i.e., the diagonal of the lower 4x4 
matrix contains the cases where the reconstruction polynomial degree matches the simulated degree. 
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Fig. 7. Influence of the dampening parameter A on the quahty of 
the cleaned maps. 



In the regime where the RMS is minimal we observe an in- 
crease of the RMS of about 2.5% (7%, 10%, 12%) for Op = 
(Op = 1, 2, 3) compared to the clean maps. Considering the fact 
that the ofl'set coefficients were drawn from a Gaussian noise 
distribution with the same standard deviation as the noise in the 
Raw data this is acceptable0. 

4. Examples 

In this section we will present two examples based on measure- 
ments with the Effelsberg 100-m telescope. 

4.1. Continuum maps (6-cm) 

The first examples is taken from to a pilot study to prove the 
feasibility of a large spectral line and continuum survey of the 
Galactic plane in the 4-8 GHz frequency range to be undertaken 
with the Jansky VLA. Supplementary Effelsberg 100-m data are 
used for short-spacings, especially for the continuum. Since cur- 
rently no receiver at the 100-m covers the 4-8 GHz range si- 
multaneously for the tests a 6-cm receiver was used providing 
500 MHz of instantane ous bandwidth. W orking in spectral-line 
mode, the new XFFTS (iKlein et al.ll2012h backend with its high 
number of spectral channels (32k) provides not only the possi- 
bility to map narrow absorption and emission line features but 
also to extract the continuum. 

For our test observations we targeted the dark cloud 
G 59.5-0.2. Two 1 x 1 deg^ fields were independently observed. 
Zig-zag on-the-ffy mapping was combined with the so-called 
position-switching technique where, after few scan lines, the 
telescope was pointed to a reference position. The latter is used 
to remove bandpass effects. To account for frequency-dependent 
system temperature, Tgys, and calibration normal, Tcai (provided 
by a n oise diode), we used the method proposed in Winkel et al. 
(|2012|) . Accordingly, every two hours a bright continuum source 
of well-known flux was observed as calibration source to deter- 
mine rcai(v). Furthermore, corrections for atmospheric opacity 
and elevation-dependent telescope eflftciency (gain curve) were 
applied. 



^ For reference: (Trms in the raw data was set to One resulting in 
c^rms ~ 0. 1 in the gridded data. 



Unfortunately, during several of the scans strong broad-band 
radio frequency interference (RFI) was entering the receiving 
system, the effect of which can be seen in Fig.[8j Here, the cal- 
ibrated data of the northern field was gridded (both coverages). 
The strong vertical artifact at / ^ 59.3° is caused by this RFI. But 
even in the remaining data there is not much to be seen, except 
for a bright continuum source in the lower left part of the figure. 

Applying the basket- weaving technique (using second-order 
polynomials) as proposed in Section [331 leads to the map shown 
in the right panel of Fig. [HI While the scan pattern is largely re- 
moved and some diff'use flux in the bottom part becomes visible, 
there is still residual RFI visible. 

To further improve the data quality, we finally masked sev- 
eral scan lines containing the strong RFI signals. This is simply 
done by removing associated rows from the solution matrix. A, 
and difference vector, D. (Of course, to calculate the correction 
map, the full matrix. A, needs to be used.) Furthermore, for the fi- 
nal gridding bad data were omitted as well. The final image com- 
posed of the northern and southern field is presented in Fig. [51 
The diffuse continuum of the dark cloud smoothly extends from 
the northern to the southern part, although no cross-calibration 
was applied between the two fields. The masking suppresses the 
residual RFI in the northern part effectively, but neglecting (a 
substantial amount of) data comes at a price: the region that was 
masked has a somewhat stronger remaining scan pattern simply 
because the solution is much less constrained there. 

4.2. Continuum maps (21 -cm, EBHIS) 

For the second ex ample we used data of the E ffelsberg-Bonn H i 
Survey ( Wink el et aLll201Ql: iKerp et al.ll201 ll) . which is a spec- 
troscopic survey of the neutral atomic hydrogen mapping the full 
northern hemisphere. It uses a seven-feed receiver with 100 MHz 
instantaneous bandwidth. A total of 14 FFT spectrometers (16k 
spectral channels each) allow one to observe the band with suflft- 
cient spectral resolution to not only detect extra-galactic sources 
out to a redshift of z ~ 0.07, but also map the Hi in the Milky 
Way. The first coverage is nearly completed. 

Here we show that the data can in principle be used to re- 
cover the continuum emission as well. For a few fields we sup- 
plemented the EBHIS data with measurements in the orthogonal 
scanning direction. Data wer e reduced with the standard EBHIS 
pipeline (IWinkel et al.ll2010h except for baseline fitting (usually 
employed to remove continuum from the spectral data). 

We extracted the continuum level at 21 -cm by calculating 
the median over the whole bandwidth. This was done to be ro- 
bust against radio frequency interference (RFI) that is present in 
roughly 10% of the channels and can greatly exceed the average 
system temperature. 

The system temperature at 21 cm is strongly dependent on 
elevation (mostly due to atmospheric opacity but with a signifi- 
cant amount of radiation from the ground). This contribution has 
to be removed first, because the basket- weaving alone will lead 
to reconstructed large-scale flux including any baseline contri- 
bution (atmosphere and ground radiation) present in both maps. 
Since the maps are observed in the equatorial coordinate system 
(scanning direction in right ascension and declination), there is 
usually a large-scale gradient present in both coverages as the 
target area moves across the (horizontal) sky. 

While the atmospheric effect can be modeled relatively eas- 
ily, the ground contribution is a more complex function of the 
telescope aperture and the shape and temperature of the hori- 
zon, we decided not to fit a certain model to the data but instead 
use a simpler approach by using a percentile filter. For each of 
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Fig. 8. Left panel: Map showing 6-cm continuum emission (after full calibration) of the northern part of the dark cloud G 59.5-0.2. 
In the left panel both coverages where simply gridded together. The right panel displays the result of our basket- weaving algorithm 
applied. Apparently a strong residual artifact remains, which we attribute to the RFI signal at / ^ 59.3° visible in the left panel. 



the seven feeds, the brightness temperatures were sorted with 
ascending elevation (Fig. [TOl upper panel) and a running per- 
centile filter, which calculates the 15% percentile from 50 val- 
ues, was applied. This filtered version, which is a lower enve- 
lope to the measured system temperatures, was then used to re- 
move the shape and amplitude of the ground level (Fig. [TOl up- 
per panel, line). Although the residual brightness temperatures 
(Fig. [TOl lower panel) have also lost information about the ab- 
solute continuum level, they should still contain information on 
the shape of the overall sky distribution, since the subtraction 
was made independently of sky position. 

After removing elevation-dependent contributions to the 
baseline level, the data were gridded; see Fig. [TT] left panel. The 
map is already quite acceptable but it can be even more improved 
with the basket- weaving method as shown in Fig.[TT]right panel. 
Especially, the scan pattern in the upper left quadrant of the im- 
age has been substantially suppressed. 

5. Summary 

We have presented a novel technique to solve the basket- weaving 
problem. Although only off'sets and not gains can be solved for 
(since it is based on linear least- squares), it has several advan- 
tages compared to previous approaches: 

- The new algorithm works on gridded data. As a main con- 
sequence it is fast (especially if strongly oversampled) and 
it can be efl&ciently applied to spectral data. Furthermore, 
arbitrary scan geometries can be wor ked with, while other 
methods rely on rectangu lar grids (e.g. JSofue & ReichI 19791: 
[Emerson & Graevd[T988h . 

- Arbitrary basis functions can be used for the off'sets in each 
scan line. In this paper we showed how to apply polynomials 
of a certain degree, but one could as well work with Fourier 
series or Legendre polynomials. These basis function can be 
arbitrarily parametrized. For example it would make sense 
to use elevation as parameter for more complicated scan ge- 
ometries since atmospheric eff'ects as one of the main con- 
tributors to scan-line artifacts in the data are mainly a func- 
tion of elevation. 

- It is very easy to handle bad data (e.g., RFI) by masking bad 
pixels in the maps and associated rows in the basket- weaving 
equation. 



- Often it is desirable to incorporate prior data, e.g., (lower- 
resolution) maps from existing surveys or from diff'erent 
wavelengths. Our technique can be easily extended to con- 
strain the solution to be consistent with such priors. 
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Fig. 9. Final 6-cm (4.8-GHz) continuum map of dark cloud 
G59.5-0.2. It was simply concatenated from the northern and 
southern fields, which were independently calibrated and cor- 
rected using the basket- weaving technique. 
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Fig. 10. Upper panel: the system temperature of a single feed 
as a function of elevation. To illustrate that multiple scan lines 
contribute at a given elevation, each scan line was color-coded. 
The line shows the result of the percentile filter (see text). Lower 
panel: the residual brightness temperatures after the subtraction 
of the percentile filter. No dependence on elevation is evident. 
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Fig. 11. Example continuum maps from the EfFelsberg-Bonn Hi Survey. Left panel: without basket- weaving. Right panel: with 
basket- weaving applied. A percentile filter was used to remove elevation-dependent contributions to the system temperature. 
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